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Abstract. - We study a Hamiltonian system describing a three spin-1/2 cluster-like interaction 
competing with an Ising-like exchange. We show that the ground state in the cluster phase 
possesses symmetry protected topological order. A continuous quantum phase transition occurs as 
result of the competition between the cluster and Ising terms. At the critical point the Hamiltonian 
is self-dual. The geometric entanglement is also studied. Our findings in one dimension corroborate 
the analysis of the two dimensional generalization of the system, indicating, at a mean field level, 
the presence of a direct transition between an antiferromagnetic and a valence bond solid ground 
state. 
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Introduction. — The Landau theory of phase transi- 
tions is at the heart of our understanding of critical phe- 
nomena [l][2] . By introducing the concept of order param- 
eter Landau's theory relates phase transitions to symme- 
try breaking. The order parameter (e.g., magnetization 
in ferromagnets) is different from zero in the less symmet- 
ric (ferromagnetic) phase while it vanishes in the para- 
magnetic phase where symmetry is restored. However, 
not all quantum phases of matter can be classified accord- 
ing to their symmetries [3] and there are quantum phase 
transitions that elude Landau's paradigm [4|[5]- These 
phases, dubbed topological phases, cannot be character- 
ized by a local order parameter. Topological phases play 
a prominent role in very diverse physical contexts, rang- 
ing from the quantum Hall effect and high Tq supercon- 
ductivity in solid state physics to confinement problems 
in QCD and string theory [3]. The continuous (direct) 
transition between a Neel antifcrromagnet and a valence 



bond solid belongs to such class of exotic phase transi- 
tions pi. The transition occurs because a new conserved 
quantity (a gauge symmetry) appears exclusively at the 
critical point, characterized by the emergence of certain 
topological defects in the zero temperature antiferromag- 
net which become energetically favorable and proliferate 
(i.e. they are deconfined). At the critical point within the 
Landau-Ginzburg scheme, such phase transitions must be 
of the first order and they cannot be continuous. 

Understanding exotic quantum phases and quantum 
phase transitions lying outside Landau's paradigm are 
among the most intriguing and controversial topics in the 
modern theory of critical phenomena [&|j8j . Exact results, 
desirable for a thorough comprehension of such phenom- 
ena, are extremely rare in the field. In this Letter we 
discuss an exotic phase transition between symmetry pro- 
tected topological order and an antiferromagnetic state 
occurring in an exactly solvable one dimensional model. 
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Based on the one dimensional theory we produce evidence 
that in two dimensions our model displays a continuous 
quantum phase transition between an antiferromagnet and 
a valence bond solid. Our model can be experimentally re- 
alized in optical lattices and ion traps (see the concluding 
remarks) and is relevant in quantum information (in cer- 
tain limits the ground state is a cluster state, central in 
one-way quantum computation I9|). 

The model. — The Hamiltonian system that we shall 
study emerges as a result of the cross fertilization between 
quantum statistical mechanics and quantum computation. 
It reads 

N N 
i=l i=l 

where erf is the Pauli matrix acting on the i-th site of a 1-D 
lattice. (Similar models were considered in 10 ll].) This 
describes the interplay between a cluster and an antifer- 
romagnetic Ising Hamiltonian. The model can be defined 
also in higher dimensions and the nature of its phase tran- 
sition is of great interest. First we focus on its properties 
in one-dimension where the model can be solved exactly. 
For A = and open boundary conditions, the ground state 
(gs) manifold is given by 



Co - span | ^ K) fe (al) 1 ]J C t |0) k, I = 0, 1 j 



where Gt = It — af 



(2) 



J ^f(^)- U ^ u ^f(^))> iS tllC 0ri " 

ented nearest neighbours of the lattice site i and cr z \0) = 
|0). Such definition is valid for a lattice of arbitrary di- 
mension. The gs manifold is thus four-fold degenerate. 
With periodic boundary conditions the gs is unique and 
is given by k = I = 0. Cluster states are a special type of 
multi-qubit graph states [9] with a computational power 
that is believed to be an important resource for measure- 
ment based quantum computation [12 1 . In this context it 
is relevant to study how the computational properties of 
the cluster states are quenched under additional pertur- 
bation. This kind of questions have been recently studied 
in a series of works [l0|[ll] where it was shown that the 
cluster phase is unstable above a critical value of certain 
one- and two-sites (Ising-like) perturbations. The compu- 
tational power of cluster states can be viewed as a quantifi- 
cation of their entanglement [l3pT 



Localizable entangle- 
ment, in particular, was calculated and it was shown that 
its range diverges 11 , a property shared with spin half 
valence bond states 15 16 . Interestingly enough the two 



states are both characterized by a hidden order of topolog- 
ical nature. Very recently it is emerging that long-range 
entanglement is a inherent property of topologically or- 
Inspired by the problem posed in [2] 



dered states 17 



we will prove that the model undergoes to a contin- 
uous transition between an Ising antiferromagnet and a 
phase characterized by the specific topological order en- 
coded in cluster states. We will show that entanglement 



captures the critical properties of this model, beyond the 
formulation of the concept of order parameter. The two 
dimensional generalization of the model will be discussed 
at the end of this letter. 

Symmetry protected topological order in cluster 
states. — Topological order is defined as a degeneracy 
of the ground state depending on the topology of the sys- 
tem. Moreover, such degeneracy must be robust under 
arbitrary local perturbations, in the sense that for a fi- 
nite system of linear size L the splitting of the degeneracy 
is exponentially small in L. Recently, it has been under- 
stood that the notion of topological order described above 
cannot be applied to ID 17|. For one dimensional sys- 
tems and therefore also for the ID cluster state, it results 
that the degeneracy is not robust under generic perturba- 
tions. Indeed, for the state ^ the four basis states are 
connected by the local operators af , cr£ and distinguished 
by the local operators crfcr% and o , J_ 1 cr| ( . Nevertheless, a 
more tenuous kind of topological order can be possible, if 
symmetries are present. If not all operators or perturba- 
tions are allowed, but only those that respect some sym- 
metry, the symmetry can protect the topological order. In 
the present case, the ground state manifold C is also de- 
fined by the algebra of the global operators (Xi,Zi) and 
( X2 , Z2) , that are defined as 
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n=0 

These operators obviously square to the identity: Xf = 
X\ = Zf = Z\ = I and, if L = 3(2fc + 1), with 
k integer, the following anticommutation relations hold: 
{Xi,Zi} — {X2,Z 2 } — and commute otherwise, thus 
realizing two copies of the Pauli algebra which constitutes 
the logical operators in the ground space Cq. The Hamil- 
tonian with A = has thus a Z2 X Z2 global symmetry 
defined by Tj = Xi + Z i: i — 1,2 that protects the topologi- 
cal degeneracy in the ground space. This means that every 
perturbation V that commutes with the symmetries is not 
able to split the degeneracy unless there is a QPT. This 
indicates that there is a whole phase in which the ground 
state degeneracy is protected. The ising term though, vi- 
olates this symmetry, so as soon as we switch on the A 
coupling the degeneracy is lifted. In the limit of large A, 
there is only the usual Z2 symmetry and the system is 
antiferromagnetic. Aspects related to the preservation of 
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topological order will be further invetigated in a forthcom- 



ing article 32 



The properties of Co reflect the emergence of low en- 



ergy states localized at the ends of an open chain 18 19 



such phenomenon can be detected by the so called string 
order parameter, originally employed to study the Hal- 
dane phase displayed by one dimensional integer-spin sys- 
tems [20]. Affleck-Kennedy-Lieb-Tasaki models provide a 
paradigmatic example in this context (see |21| for a recent 
reference), with symmetry protected topological ordered 
ground states 17,21,22. We observe that (p| is sym- 
metric for alternating spin flips, J2k(~ (^o\^1\Cq) = 0, 
and any local order is absent. The string order parameter 
O z = i-) N - 2 (<rlUjJi°j°N) # 0. reflects the symme- 
try breaking of the state by n rotations about the y and 
the z axes. It can be shown that O z = 1 for the cluster 
state, i.e. the gs of |l| at A = 0. The gs of the anti- 
ferromagnetic Ising Hamiltonian \Iq) corresponds to the 
A — > 00 limit, instead. The symmetry of the system is 
spontaneously broken, so that the staggered magnetiza- 
tion M = J2k(~ l) fe (^ol cr fc \Io) is an extensive quantity in 
the thermodynamic limit. 

Duality mapping. — Duality transformations pro- 
vide a powerful tool to extract properties of the phase 
diagram beyond the perturbative regions 
apply the duality transformation 



23 . We now 
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n 



(4) 



to our Hamiltonian. Observe that Eq. Q is a unitary 
transformation. 

We first consider the case A = in Eq. ([lj. For open 
boundary condition, (Tq = crj^+i = 1, the cluster Hamilto- 
nian if (A = 0) is transformed into an Ising Hamiltonian 
with a boundary term: 



if (A = 0) if 7 



N-l 

»=i 



B, 



(5) 



where B — i(J\i=i M* )Mi ■ Being unitary, the dual map- 
ping Q factors out the correlation terms into a bound- 
ary term. The gs of ^ is a superposition of states with 
anti-ferromagnetic order, \G) = (\Iq) + \Io))/y2, with 

Uo) = |+»,-»,+y,-V-), l^o) = B\I ), ^\±y) = ±\±y), 

and Hj\G) = —N\G), where — N is the non-degenerate 
lowest eigenvalue of the Hamiltonian. \G) can be trans- 
formed into a GHZ state |GHZ) = (|+) 8JV + 1-)®^)/^ 
by local unitary transformations. Therefore, after the 
dual transformation, the cluster state is transformed into a 
GHZ state and the unit of entanglement is degraded from 
N/2 to 1 (Til . On the other hand, for periodic boundary 
conditions, Eq. Q yields 



if (A = 0) if 



N-2 
i=l 



B 



(6) 



where B p = (H i=1 ^f)(i^ z N - I^n-iVn) ■ The § s is a S ain 
a GHZ state \G p ) = (\I )n-i + |Iq}jv-i) ® \+x)/V2 and 
the same considerations apply. 

We now tackle the case A 7^ 0. The dual model of |l]) is 

N-2 N-l 

if (A) dual = Yl - A £ liU^H-i - B, (7) 

i=l i=l 



where /ijj = 1 and B = (]jf =1 (if ) (ifif^N ~ Vn-i^n + 
ifJ-N-i^N^i)- I n the thermodynamic limit N — > 00, the 
duality relation reads if(A)d U ai = Aif(A _1 ), since the 
boundary can be neglected [24]. The self-duality relation 
at A = 1 signals the critical point. 

By resorting to the duality properties ^ the local or- 
der parameter of the antiferromagnet is mapped into the 
string order parameter of the cluster states (see also [Tol 



11 ). In particular we note that the Si( cr f-i cr f cr f+i)/-^ 
cannot be taken as a valid order parameter because it en- 
joys the same symmetries of the Hamiltonian. 

Fermion representation and QPT. The model 
defined here can be diagonalized by resorting to a Jordan- 



Wigner transformation c\ = YimLi 
(spinless fermions), yielding 



. CT m CT z + > c \ c i 



1/2 



N 



H(X) = X)( c i-i-Ci-i)(4fi + Q+i) 

2 = 1 

JV 

+ A 5>l +ci)(c\ +1 -c l+1 ), (8) 

Z=l 

with boundary conditions c/v+z = TQ- N is assumed 
to be even. We note that the three spin interaction of 
the cluster Hamiltonian is reflected in the next-nearest 
neighbor interaction in the fermionized Hamiltonian 
The dispersion is obtained by Fourier transformation c\ = 
2^ fe e^ lkl l N bkj \f~N followed by a Bogoliubov transforma- 
tion b k = u k j k + v k T_ k , yielding: 

if (A) = £A fe ( 7j U-l/2), 



A fe = yj{\ 2 + 1) - 2Acos (Qirk/N), 



(9) 



with u k = z£, v k = -isiga(S k )z k , = y/(T± e k /A k )/2, 
e k = [cos(4wk/N) - A cos(27rfc/JV)], 6 k = [sin(47rfc/iV) + 
Asin(27rfc/A)] and -(N/2 + l)/2 < k < -1, N = 4n + 2 
[-(N - l)/2 < k < -1/2, N = An] for [anti] periodic 
boundary conditions. The gs is identified as the vacuum 
of the Bogoliubov operators 7fc|fi) = 0, Vfc. In terms of 
the b k 's it can be written as the BCS ground state 



\V)=I[(uk+V k blbl k ) \vac). 



(10) 



From simple inspection of the dispersion law of Afc, one 
infers that the system displays a second order quantum 
phase transition at A = 1 with critical indices z = v = 1. 
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Notice also that even though any two-spin correlations 
vanish in the cluster states [l2|[25] , the fermionic pair cor- 
relations pattern is not trivial. 

Entanglement is a key quantity in the study of many- 
body systems 26 . We now focus on how maximal entan- 



glement in the cluster state 13 is demoted by nearest- 
neighbor interaction, as A increases in model ( Tj). We 
consider a geometric measure of entanglement 27 that 
provides a global characterization of the entanglement in 
quantum many-body systems 



e (|^)) = -log 2 



:|<S(W)| 2 



(11) 



where \S(9)) is the closest separable state. iV-partite pure 
separable states are characterized by 2N real parameters, 

N 

\S({9A)) = II ( cos ^ + ^ 8in VJ) I t> 0JV - (12) 
J'=l 

that can be restricted to four because of the dimerized 
structure of the eigenstate of the Hamiltonian. 
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entanglement is very small for A > 1, it does not vanish. 
Observe that the convexity of the geometric entanglement 
changes from positive to negative near the critical point. 
As N is increased, the point of convexity change moves to- 
ward A = 1. This is scrutinized in the left panel, where one 
sees that the phase transition is detected by the derivative 
of the geometric entanglement, that diverges at A = 1. 

d-dimensional systems. — The model defined by the 
Hamiltonian ([I]) can be generalized to d-dimensions, the 
cluster term involving the product of 2d + 1 spin operators 
(we consider for simplicity a hyper-cubic lattice). We do 
not rely on exact results in this case; nevertheless we will 
show that the system has a critical value of the parameter 
A below which the antiferromagnetic order vanishes. This 
argument, together with the remarkable observation due 
to Cirac and Verstraete [29] that cluster states can be real- 
ized as valence bond states, leads us to conclude that our 
model displays a direct transition between an antiferro- 
magnet and a valence bond state, for d > 1. We comment 
that the treatment we are going to employ is mean-field 
in nature and although it is not expected to capture the 
correct critical behavior, it suffices to establish the very 
existence of the transition point. 

The mean-field analysis [3l] proceeds after a rotation 
of the Hamiltonian to map the antiferromagnetic to ferro- 
magnetic Ising couplings. The approximated Hamiltonian 
reads: 



MF 



n < 



2d\iP a i 



V (13) 



where tp = (of) is the (real and spatially uniform) or- 
der parameter and 2d the coordination number of the lat- 
tice. By assuming that the order parameter vanishes at 
the critical point, the self-consistency equation leads to 
the condition 



1 = 4dA. 



E 



dTe-^-^\(C\af\C m }\ 2 , (14) 



Fig. 1: Left panel: The first derivative of the geometric entan- 
glement captures the critical point at A = 1; N ranges from 
32 (black-dashed line) to 1000 (turquoise-triangle line). Right 
upper pannel: Normalized geometric gs entanglement e/N as 
function of A. The decay of entanglement in the Ising phase 
is shown in the lower right panel: squares and crosses are the 
data for e/N at the critical point and in the antiferromagnetic 
phase, respectively (an artificial offset has been applied for a 
convenient display). 



The geometric entanglement ( 11 ) is plotted in Fig. [I] In 
the right upper panel the normalized geometrical entan- 
glement is plotted as a function of A for different values of 
N. For A = 0, the geometric entanglement is always N/2. 
In the cluster phase (A < 1) entanglement is very weakly 
dependent on the size of the system, while in the antifer- 
romagnetic phase (A > 1) it decreases with N. Moreover, 
as shown in the right lower panel, although the normalized 



where \C m ) and e m are the m-th eigenstate and the energy 
level of the cluster Hamiltonian. A direct evaluation leads 
to A c = 1 in any dimension, marking the quench of the 
magnetic order. The independence of the critical value on 
the dimensionality reflects the peculiar multi-spin interac- 
tion of the cluster couplings. 

Conclusions. — Wc studied the phase diagram and 
the geometric entanglement of the (spin 1/2) cluster 
Hamiltonian in the presence of an additional Ising inter- 
action. We showed that Cluster states have symmetry 
protected topological order, which can be quenched to a 
long range order by varying the two spin interaction in a 
continuous way. The one dimensional model we studied in 
this paper is equivalent to a free fermion Hamiltonian, and 
therefore it can be solved exactly. The Hamiltonian enjoys 
a duality symmetry, mapping the three point correlations 
into nearest neighbor interactions; as a consequence the 
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entanglement in the cluster state results to be dual to the 
entanglement encoded into the Ising gs state. The system 
is self-dual at the critical point. We have shown how the 
geometric entanglement characterizes the two phases and 
the QPT. The exact solution in one-dimension together 
with the argument we produced above for higher dimen- 
sions provides, in our opinion, a convincing scenario for 
the transition from Neel order to a cluster state order- 
ing of topological nature. The critical properties of the 
transition could be investigated by exploiting the scheme 
developed in 
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These results are also relevant for quantum informa- 
tion technology. It would be valuable, for example, to ex- 
ploit the Ising-cluster state duality relation as a resource 
for computational power in quantum algorithms. Finally, 
we observe that the cluster Hamiltonian can be realized 
in a triangular optical lattice of two atomic species 



2cS 



The three-spin interactions arises from a kind of density- 
dependent tunneling of the atomic species. Our results 
can therefore be tested in optical lattices and in trapped 
ions, where this model can be experimentally realized. 
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